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Abstract 

Coronaviruses are single-stranded positive-sense RNA viruses associated with important avian diseases. Their relatively high rates of 
mutation and recombination frequencies allow them to adapt to new hosts and ecological niches. Although Brazil has 1896 of global 
avian species diversity, studies regarding the presence of avian viral diseases in wild birds in South America are scarce. In this study, 
we performed a retrospective analysis of the presence of CoVs in 746 wild birds. Oropharyngeal and cloacal swabs were obtained and 
placed together in vials containing VTM transport medium collected in different regions of Brazil between 2006 and 2013. Screening 
for viral nucleic acid was performed using conventional RT-PCR and pancoronavirus nested PCR. Positive samples were character- 
ized by partial sequencing of the RNA-dependent RNA polymerase (RdRp) gene, and ensuing phylogenetic analysis was performed 
to investigate the association between virus epidemiology and bird migration routes. Coronavirus RNA were detected and sequenced 
from six samples, in which three were related to gammacoronaviruses group and the other three to deltacoronavirus group. Our study 
documents the presence of CoVs related to avian gamma- and deltacoronaviruses circulating in both urban- and poultry-farm regions 
of Brazil, implicating wild birds as potential carriers of CoVs which may represent a risk to poultry farms and public health in Brazil. 
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Introduction Infectious bronchitis virus (IBV), known as coronavirus of 

chickens, causes important economic losses to the poultry in- 
Coronaviruses are RNA single-stranded positive polarity virus dustry. Studies about the presence of CoVs in wild birds, es- 
associated with important disease in avian poultry [1]. pecially IBV, in South America are scarce [2]. Coronaviruses 
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(CoVs) possess the largest genomes among RNA viruses and 
have relatively high mutation rates and recombination frequen- 
cies when compared with other single-strand RNA viruses, 
which may allow them to adapt to new hosts and ecological 
niches [3—5]. These characteristics are exemplified by the two 
emerged highly pathogenic betacoronaviruses, CoV-SARS 
(severe acute respiratory syndrome) and the CoV-MERS 
(Middle East respiratory syndrome) [6]. 

Even though avian coronaviruses with zoonotic potential have 
not been reported yet, many species of wild and domestic birds 
can carry CoVs responsible for substantial economic loss within 
the poultry industry such as the gammacoronavirus infectious 
bronchitis virus (IBV). This virus affects the gastrointestinal, 
respiratory, urinary, and reproductive systems and may lead to 
mortality rates as high as 30% and sudden drops in egg produc- 
tion [7]. Moreover, many novel CoVs genetically distinct from 
IBV have been recently identified in different families of wild 
and domestic birds and mammals but their potential economic 
losses remain unknown. Recently, Duráes-Carvalho et al. (2015), 
described the detection of delta- and betacoronaviruses in captive 
wild birds, from a rehabilitation center, originally captured in São 
Paulo and Mato Grosso do Sul state with close evolutionary 
relationships with CoVs isolated from mammals [8]. 

Wild migratory birds may facilitate the spread of viruses, as 
demonstrated by the surveillance of avian influenza virus [9, 
10] and West Nile virus [11—13], but the circulation of CoV in 
avian migratory routes in South America, more specifically in 
Brazil, still need to be elucidated. The potential diversity of 
CoV in Brazilian bird fauna as consequence of the bird species 
diversity (18%), together with the large number of endangered 
species and the high intensity of poultry farms in Brazil, indi- 
cates that birds from Brazilian biomes are important subjects 
for CoV surveillance [2, 14]. 

In this study, we performed a retrospective analysis, mo- 
lecular detection, and phylogenetic analysis of avian 
coronaviruses using oropharyngeal and cloacal swab samples 
from wild and domestic birds collected between 2006 and 
2013 in different regions of Brazil to investigate the associa- 
tion between virus epidemiology, bird migration routes, and 
the proximity of urban areas and poultry farms. 


Materials and methods 
Sample collection 


In total, 746 samples (Table 1) of birds captured between 2006 
and 2013 in five expeditions to north, northeast, midwest, south, 
and southeast Brazil (Fig. 1) were analyzed. For selection, the 
migratory routes were investigated through a map analysis and 
experience of ornithologists [15]. The samples were collected by 
field team of the Laboratory Clinical and Molecular Virology, 
University of São Paulo (USP), Sao Paulo, Brazil, with support 
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of the Wildlife Department of São Paulo (DEPAVE); University 
of Vale do Rio dos Sinos (UNISINOS); University of 
Pernambuco (UFRPE); National Center for Research and 
Conservation of Wild Birds (CEMAVE); Laboratory of Bird 
Ecology, Institute of Biosciences (UFMT); and Laboratory of 
Virology and Rickettsioses, Veterinary Hospital (UFMT). 
Samples were collected from wild birds after capture with mist 
nets, strategically placed between 30 and 250 cm of the ground, 
and by hand nets on an open field in public parks. The field 
works were made according to the permission of the Brazilian 
Institute of Environment and Renewable Natural Resources 
(IBAMA) and National Center for Research and Conservation 
of Wild Birds (CEMAVE) developed and followed protocols 
approved by the Ethics Committee on Animal Experimentation 
of the Institute of Biomedical Sciences, University of Sao Paulo. 
Sampling from oropharyngeal and cloacal swabs was obtained 
and placed together in vials containing VTM transport medium 
[16]. All birds were then released after sampling. 


RNA extraction and reverse transcription 


Viral RNA was extracted using the MagMax TM-96 RNA 
Isolation Kit (Ambion, Austin, TX, USA) according to the 
manufacturer’s instructions. The extracted RNA was eluted 
in 80 uL of RNase-free water and stored at — 80 °C until 
processed. The cDNA was transcribed using the High- 
Capacity cDNA Archive Kit (Applied Biosystems, Foster 
City, CA, USA) following the manufacturer's instructions. 
The RT reaction was performed by adding, to each reaction, 
10 uL of total RNA, 2 uL of RT Buffer 10X, 2 uL of Random 
Primers (50 ng/uL) (Thermo Fisher Scientific, Waltham, MA, 
USA), 4.2 uL of DEPC water, and 50 U of MultiScribeTM 
Reverse Transcriptase (Thermo Fisher Scientific). The RT step 
conditions were 25 °C for 10 min, 37 °C for 120 min, and 
85 °C for 5 min, and cooling to 4 °C. 


Coronavirus detection 


Polymerase chain reaction (PCR) followed by nested PCR 
were performed according to the protocol for the detection 
of avian CoVs developed by Chu et al. (2011) [17], which 
targets the RNA-dependent RNA polymerase (RdRp) gene 
in the ORFI area of the viral genome. The resulting approx- 
imately 600-bp PCR fragment and 440-bp nested PCR frag- 
ment were separated by agarose gel [1.5%] electrophoresis 
and viewed under UV light after staining with ethidium bro- 
mide (25 ug/100 mL). 


Purification and sequencing 
Samples were purified using the ExoSap-IT enzyme (GE 


Healthcare, Chalfont St. Giles, UK) and GeneJET™ Gel 
Extraction Kit (Thermo Fisher Scientific) according to the 
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Table 1 Surveillance of Coronavirus in resident and wild birds in Brazil by molecular tests. Bold represents positive samples 


Site Date Scientific name Sample* No. pos./No. birds 
Ilha de Marajó, PA Jul. 2006 Cairina mochata C and O 0/42 
Gallus gallus C and O 0/36 
Meliagres gallopavo C and O 0/15 
Ajuruteua, PA Nov.2006 Actitis macularius C and O 0/5 
Calidris pusilla C and O 0/2 
Charadrius semipalmatus C and O 0/1 
Dendrocygna autumnalis C and O 0/18 
Arenaria interpres C and O 0/9 
Dendrocygna viduata C and O 0/56 
Ilha de Canelas, PA Nov. 2006 Actitis macularius C and O 0/10 
Calidris canutus C and O 0/1 
Calidris pusilla C and O 0/3 
Calidris minutilla C and O 0/1 
Sterna hirundo C and O 0/4 
Larus sp. C and O 0/1 
Nov. 2008 Actitis macularius C and O 0/31 
Arenaria interpres C and O 0/22 
Calidris alba C and O 0/2 
Limnodromus griseus C and O 0/9 
Tringa melanoleuca C and O 0/1 
Charadrius semipalmatus C and O 0/1 
Gelochelidon nilotica C and O 0/1 
Bahia de São Marco, MA Apr. 2007 Dendrocygna viduata C and O 0/27 
Netta erythrophthalma C and O 0/2 
Dendrocygna autumnalis C and O 0/48 
Gallinula chloropus C and O 0/2 
Amazonetta brasiliensis C and O 0/14 
Anas bahamensis C and O 0/13 
Meliagres gallopavo C and O 0/1 
Lagoa do Peixe National Park, RS Nov. 2009 Rynchops niger C and O 1/63 
Calidris fuscicolis C and O 1/63 
Haematopus palliatus C and O 0/6 
Himantopus melanurus C and O 0/2 
Pluvialis dominica C and O 0/7 
Arenaria interpres C and O 0/4 
Calidris canutus C and O 0/3 
Calidris alba C and O 1/6 
Calidris pusilla C and O 0/1 
Nycticriphes semicollaris C and O 0/1 
Fulica leucoptera C and O 0/1 
Puffinus gravis C and O 0/1 
Puffinus puffinus C and O 0/3 
Tringa flavipes C and O 0/2 
Sterna hirundo C and O 0/15 
Phalacrocorax brasilianus C and O 0/2 
Vanellus chilensis C and O 0/1 
Elaenia mesoleuca C and O 0/1 
Pluvialis squatarola C and O 0/1 
Tryngites subruficollis C and O 0/7 
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Table 1 (continued) 
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Site Date Scientific name Sample* No. pos./No. birds 
Calidris melanotus C and O 0/1 
Sterna vitatta C and O 0/1 
Sao José do Egito Caatinga, PE Oct. 2010 Thaminophilus capstratus C and O 0/8 
Hylophilus amaurocephalus C and O 0/10 
Sarkesphorus cristapus C and O 0/8 
Hemitricus margaritaseiventer C and O 0/4 
Coryphospingus pileatus C and O 0/10 
Anopetia gounellei C and O 0/1 
Casiornis fuscus C and O 0/1 
Gyolophylox hellmayri C and O 0/3 
Tolmyias flaviventris C and O 0/2 
Heliomaster squamosus C and O 0/2 
Cychlars gujanensis C and O 0/1 
Cyanoloxia brissoni C and O 0/1 
Formicivora melanogastrer C and O 0/8 
Cantorchilus longirostris C and O 0/1 
Chlorostilbon lucidus C and O 0/1 
Coereba flaveola C and O 0/1 
Polioptila plumbea C and O 0/1 
Synalax frontalis C and O 0/1 
Pantanal Poconé, MT Nov. 2012 Ramphocelus carbo C and O 0/10 
Hemitriccus margaritaceiventer C and O 0/6 
Cnemotriccus fuscatus C and O 0/11 
Paroaria capitata C and O 0/2 
Picumnus albosquamatus C and O 0/3 
Gallus gallus C and O 0/9 
Furnarius rufus C and O 0/8 
Myiarchus swainsonii C and O 0/2 
Ibirapuera Park, SP Jun. 2013 Anser cygnoides C and O 3/81 
Total 63 species 6/746 


*Sample type: C = cloacal swabs and O = oral swabs 


manufacturer’s instructions. The Sanger sequencing reaction 
was prepared using Big Dye Terminator v3.1 Cycle 
Sequencing Kit (Life Technologies, Foster City, CA, USA), 
followed by purification using BigDye X Terminator® 
Purification Kit (Applied Biosystems), and analysis in an au- 
tomated ABI PRISM 3130XL DNA Sequencer (Applied 
Biosystems). 


Genetic characterization and phylogenetic tree 
of the RdRp gene 


Generated chromatograms containing nucleotide sequences 
were visualized using Bioedit Sequence Alignment Editor 
(v.7.1.5.0). Nucleotide sequences were aligned using CLC 
sequencer WorkBench version 7.6 (QIAGEN Aarhus, 
Denmark), and the nucleotide percentage similarities were 
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calculated using the MegAlign program of the DNASTAR 
package (DNASTAR, Inc., Madison, WI, USA). The similar- 
ities with existing sequences were determined using the Basic 
Local Alignment Search Tool (BLAST). Phylogenetic analy- 
sis was performed using the neighbor-joining method (NJ) 
implemented in the MEGA 7.0.26 [18]. The nucleotide sub- 
stitution model used for phylogenetic reconstruction was the 
p-distance model. The bootstrap values > 55% were obtained 
in the analysis of 10,000 replicates and are presented at the 
branching points. 


Fig. 1 Location of sampling sites with positive results for avian CoVs > 
and their relationship with the main bird migratory routes and 
congregation sites in North and South America. Congregation sites 
include reproduction, wintering, and tagging areas of boreal migrant 
birds (modified from http://www. icmbio.gov.br/cemave/) 


Braz J Microbiol 


A Delaware Bay 


JIha de Marajó- PA 
Parasinho Island-PA Wha de Canelas- PA 
Maranháo 
Récesses -MA 


Sho José 
dofgito- PE 


Pantanal-MT 


Legend riesen 


mm Migratory Route 
Reproduction Site 
—— Wintering Site 
BE Staging Site 
@ Sampling Site With Positive Results 


San Antonio Oeste 


© Sampling Site Without Positive Results 


A Springer 


Table 2 Sampling location, total 
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No. of samples Positive samples Year 
93 - 2006 
91 - 2006 
67 — 2006—2008 
107 = 2007 
192 3 2009 
64 - 2010 
51 - 2012 
81 3 2013 
746 6 2006-2013 


FA Para, SP Sao Paulo, MT Mato Grosso, BA Bahia, RS Rio Grande do Sul 


samples, number of positive Location 

samples obtained from birds 

collected across Brazil from 2006 Ilha de Marajó/PA 

to 2013 Ajuruteua/PA 
Ilha de Canelas/PA 
Bahia de Sáo Marco/MA 
Lagoa do Peixe National Park/RS 
Sao José do Egito/PE 
Pantanal Poconé/MT 
Ibirapuera Park/SP 
Total 

Results 


The CoV screening protocol, based in a nested PCR assay 
with degenerate consensus primers targeting the conserved 
region of RNA-dependent RNA polymerase (RdRp) gene, 
was selected for its greater sensitivity and CoV genera ampli- 
tude detection. Ofthe 746 bird samples collected from 2006 to 
2013 in different regions of Brazil and analyzed by conven- 
tional RT-PCR, six samples (0.896 of total) from two different 
regions were confirmed by DNA sequencing (Table 2). Our 
study describes the detection of samples related to avian 
gamma- and deltacoronaviruses in birds collected in southeast 
and southern region of Brazil, more specifically in São Paulo 
and Rio Grande do Sul state (Table 3). The presence of CoV 
RNA was detected in three Chinese geese (Anser cygnoides) 
sampled in Ibirapuera Park, downtown Sao Paulo, southeast 
Brazil. One sample has shown similarities (93%) with se- 
quences from viruses in the gammacoronavirus group, which 
are related to migratory birds only, whereas the other two 
samples (84%) belong to the deltacoronavirus group. Even 
though these sequences are closely related to those of wild 


Table 3 
collected across Brazil from 2006 to 2013 


Sample Species Location Group 

DPV 5 Anser cygnoides Ibirapuera Park/SP — Deltacoronavirus 

DPV 10 Anser cygnoides Ibirapuera Park/SP — Deltacoronavirus 

DPV 16 Anser cygnoides Ibirapuera Park/SP = Gammacoronavirus 

PNLPIOO Calidris alba Lagoa do Peixe Gammacoronavirus 
National Park/RS 

PNLPII5 Rynchops niger Lagoa do Peixe Deltacoronavirus 
National Park/RS 

PNLP159 Calidris fuscicollis Lagoa do Peixe Gammacoronavirus 
National Park/RS 


SP Sao Paulo, RS Rio Grande do Sul 
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bird CoVs, this group also includes ferret, pig, and leopard 
CoVs, demonstrating the potential for spillover events be- 
tween coronaviruses from bird and mammal hosts. In the same 
way, in southern Brazil, three samples positive for CoVs were 
obtained from migratory birds captured in wintering areas at 
Lagoa do Peixe National Park, state of Rio Grande do Sul. 
Phylogenetic analysis showed that these samples (PNLP 100, 
PNLP 115, and PNLP 159) belonged to two CoV genera. 
Samples related to gammacoronaviruses were obtained from 
one Calidris alba and one Calidris fuscicollis individual, re- 
spectively, related a restricted to sandpiper species (Fig. 2). 
Samples related to deltacoronavirus were obtained from a 
black skimmer (Rynchops niger) carrying a strain that clus- 
tered, in phylogenetic analysis, closer to the clade of viruses 
found in Ibirapuera Park (Fig. 3). 


Discussion 


The nested PCR assay with degenerate consensus primers was 
used to many studies. This assay has been used for the 


Sample identification, bird species, place of capture, and coronavirus group of avian coronavirus sequences obtained from bird samples 


Accession 
number 


Year (GenBank reference) Similarity (96) 


2013 KU321643 Wigeon coronavirus HKU20 84 
strain HKU20-9243 (2012) 

Wigeon coronavirus HKU20 84 
strain HKU20-9243 (2012) 

Avian coronavirus isolate 93 


7 M/Porphyrula alleni (2011) 


2013 KU321644 


2013 KU321645 


2009 KU321640 Rock Sandpiper Coronavirus 99 
(2010) 

2009 KU321641 Chroicocephalus 81 
ridibundus/Finland (2013) 

2009 KU321642 Western Sandpiper 99 


Coronavirus/KR-28 (2010) 
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100  GQ427173.1/Turkey/coronavirus/TCoV/VA-74/03 


i E L| EUO95850.1/Turkey/coronavirus/MG 10 
E L— KJ690953.1/Mallard/coronavirus/ACoV/Mallard/PL-MW345/2009 
C KJ135013.1/Chicken/infectious bronchitis virus/IBVUkr27-11 


PC 


I Lagoa do Peixe National Park 
Wi Ibirapuera Park 


]$9 


Bi Positive Control- IBV 


Gammacoronavirus 


| 100 


D 

= 109 Positive Control- Avian Infectious bronchitis virus strain Ma5 
KJ425488.1/Chicken/infectious_bronchitis_virus/ck/CH/LHB/110825 

KJ690955. 1/Bean_Goose/coronavirus/ACoV/Bean_goose/PL-MW435/2009 


KM093893.1/Porphyrula_alleni/coronavirus/CV 11/2011 

it KM093897.1/Bulbucus_ibis/coronavirus/6cvPF 12/2011 

o KM093895.1/Gallinago macrodactyla/coronavirus/AcvPF 7/2011 

KM093899.1/Bulbucus ibis/coronavirus/9cvPF 17/2011 

i KM093879.1/Charadrius pecuarius/coronavirus/S22/2011 

| KM093885.1/Gallinula_chloropus/coronavirus/4P/2011 

KM093883.1/Charadrius_pecuarius/coronavirus/2P/2011 

[Hi KM093882.1/Rostratula benghalensis/coronavirus/1P/2011 

KM093881.1/Gallinula chloropus/coronavirus/S23/2011 

! KM093880.1/Anas_hottentota/coronavirus/S21//2011 
PNLP100/Calidris_alba/coronavirus/Brasil/2009 

95] PNLP159/Calidris_fuscicollis/coronavirus/Brasil/2009 


ja |_|! GU396675.1/Western_sandpiper/coronavirus/KR-28 


Ju dE 
] GU396685.1/Rock sandpiper/coronavirus/CIR-665821 


zm KM093892. 1/Porphyrula alleni/coronavirus/7M/Porphyrula alleni/2011 


100 


Deltacoronavirus 


CE 
— 
0.06 


Fig. 2 Phylogenetic analysis of partial RdRp gene sequence of avian 
gammacoronaviruses. Phylogenetic tree of partial sequences of the 
RdRp-encoding region constructed using the neighbor-joining method. 
In green, CoVs related to gammacoronaviruses sequences detected in 
samples from Lagoa do Peixe National Park, state of Rio Grande do 


identification and differentiation of CoV in the four genera 
described to date (alpha-, beta-, gamma-, and 
deltacoronavirus) [17, 19, 20]. Similarly to our report, most 
of avian CoV diversity studies detected the same coronavirus 
genera [2, 17, 21, 22], with the exception of a report of 
betacoronavirus detection in distinct birds species from a re- 
habilitation center in Brazil [8]. Further studies are needed to 
understand the circulation of betacoronavirus in birds. 
Ibirapuera is the main park in Sao Paulo, receiving on 
average 150,000 visitors in a single weekend. Wild birds 
make up a large portion of the species of animals that in- 
habit the park (ca. 142 species). The birds are kept free and 
cohabiting with the human population and other wild species 
(resident or visitors). There are another 21 species of wild 
animals, including bats and skunks, and the park also allows 
pets [23]. Such close contact between so many different 
species may cause interspecific transmission and 


KM093890.1/Porphyrula_alleni/coronavirus/8CV/Porphyrula_alleni/2011 
100| KM093901.1/Bulbucus_ibis/coronavirus/6cvPF 11/Bubulcusibis/2011 
KM093891.1/Porphyrula alleni/coronavirus/5P/Porphyrula alleni/2011 
DPV 16/Anser cygnoides/coronavirus/Brasil/2013 
ss [— GU396690.1/Snow  goose/coronavirus/WIR-159 
[ousonerr uen nosotras 
GU396676.1/Brent  goose/coronavirus/KR-70 
L- JN788879.1/Dendrocygna javanica/coronavirus/KH08-0852/Dendrocygna 


Sul, southern Brazil. In red, CoVs related to gammacoronaviruses 
sequences detected in samples from Ibirapuera Park, Sào Paulo, southeast 
Brazil. Bootstrap values > 5596 were obtained in the analysis of 10,000 
replicates and are presented at the branching points 


recombination events, leading to the appearance of emerging 
viruses [24]. The three goose samples were collected from 
geese that live in this same environment in contact with wild 
birds that use the park as a stopping point. On the other 
hand, three migratory birds captured in wintering areas at 
Lagoa do Peixe National Park-RS have shown positive re- 
sults of CoVs too. Sandpipers are boreal migrants that mi- 
grate south from North America in the winter [25]. A CoV 
related to deltacoronavirus group was obtained from a black 
skimmer (Rynchops niger), thus highlighting the possibility 
that migratory birds may be carrying these viruses across 
migratory routes in Brazil [9]. The Lagoa do Peixe region 
is an important migratory route [25], and it is located close 
(400 miles) from one of the most important region of poultry 
farm production of Brazil, more specifically region of Porto 
Alegre in Rio Grande do Sul State [26]. The proximity be- 
tween the migratory routes of wild birds infected by 
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Fig. 3 Phylogenetic analysis of partial RdRp gene sequence of avian 
deltacoronaviruses. Phylogenetic tree of partial sequences of the RdRp- 
encoding region constructed using the neighbor-joining method. In green, 
CoVs related to deltacoronavirus sequences detected in samples from 
Lagoa do Peixe National Park, state of Rio Grande do Sul, southern 


gammacoronavirus shows the importance of a continuous 
surveillance and potential contact control between poultry 
and migratory birds as measures to avoid any chance of viral 
transmission between these animals. 

The CoVs related to deltacoronaviruses group (DPV 5 and 
DPV 10) and gammacoronavirus group (DPV 16) from 
Ibirapuera Park showed the highest similarity to a CoV found 
by Lau et al. (2012) [27] in Guangdong province (84%), 
China, and by Lima et al. (2015) [28] (93% nt identity) in a 
migratory bird (Porphyrula alleni) from Madagascar in 2011, 
respectively (Table 2). The gammacoronaviruses found at 
Lagoa do Peixe National Park were more similar to those 
found by Muradrasoli et al. (2010) [29] in sandpipers 
(Calidris spp.) from the Bering Strait (99% nt identity). 
However, the deltacoronavirus was more similar to the strain 
isolated from a black-headed gull (Chroicocephalus 
ridibundus) by Hepojoki et al. (2017) in Finland (95% nt 
identity) (Table 3). Interestingly, this data were reported about 
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a KJ601779.1/Swine/Deltacoronavirus/PDCoV/USA/Illinois136/2014 


KJ481931.1/Swine/Deltacoronavirus/PDCoV/USA/Illinois121/2014 
100 | KP981395.1/Porcine/Deltacoronavirus/USA/IL/2014/026PDV 
KJ769231.1/Swine/Deltacoronavirus/OhioCVM1/2014 
KJ601780.1/Swine/Deltacoronavirus/PDCoV/USA/Ohi0137/2014 
KJ567050.1/Porcine/Deltacoronavirus/8734/USA-IA/2014 
KMO012168.1/Porcine/coronavirus/HKU 15 
EF584908.1/Leopard cat/coronavirus/Guangxi/F230/2006 
KP757890.1/Porcine/Deltacoronavirus/CHN-AH-2004 
10 | EF584911.1/Chinese_ferret_badger/coronavirus/Guangxi/F247/2006 
JQ065042.1/Porcine/coronavirus/HKU 15 
81 KJ584356.1/Porcine/coronavirus/HKU15 
100 | KM820765.1/Porcine/Deltacoronavirus/KNU 14-04 
JQ065045.1/Sparrow/coronavirus/HKU17 
100 ES 
FJ376620.1/Bulbul/coronavirus/HKU1 1-796 
FJ376621.1/Thrush/coronavirus/HKU 12-600 
58 JQ065044.1/White-eye/coronavirus/HKU 16 
JX548304.1/Shorebird/coronavirus/85 
KX588674.1/Chroicocephalus ridibundus/Finland/2013 
PNLP115/Rynchops niger/coronavirus/Brasil 
JQ065049.1/Common-moorhen/coronavirus/HKU21 
FJ376622.1/Munia/coronavirus/HKU 13-3514 
JQ065046.1/Magpie-robin/coronavirus/HKU18 
100 | DPV. 5/Anser cygnoides/coronavirus/Brasil 
65 DPV 10/Anser cygnoides/coronavirus/Brasil 
JQ065048.1/Wigeon/coronavirus/HKU20 


100 


Brazil. In red, CoVs related to deltacoronavirus sequences detected in 
samples from Ibirapuera Park, São Paulo, southeast Brazil. Bootstrap 
values > 5596 were obtained in the analysis of 10,000 replicates and are 
presented at the branching points 


sampling collected during the same years (2010-2013) by the 
same way, using RT-PCR targeting the RdRp gene to detect 
the presence of CoV in wild birds from Finland. The preva- 
lence was found around 5.496 (51/939), which corroborates 
with our findings. Hepojoki et al. (2017) have shown positive 
samples to diverse gamma- and deltacoronavirus too [22]. 

In South America, there is a lack of information about 
coronavirus in wild birds. These data are interesting, once 
that, present birds infected in central park in Sao Paulo 
city and wild birds through migratory routes. In Hong 
Kong, Chu et al. (2011) detected a high prevalence 
(12.596) in aquatic wild birds, reporting gamma and 
deltacoronavirus, being the first one more prevalent [17]. 
Similar studies screened before in Sweden (2016), with 
the number of sampling similar to our study, and reported 
a high prevalence (18.7%) in wild birds, with emphasis 
for Anseriformes that had the highest prevalence (51%), 
all belonging to the gammacoronavirus group [30]. On the 
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other hand, Jordan et al. (2015) analyzed 303 samples 
from cloacal swabs, using universal CoV primers, and 
found only one duck coronavirus, belonging to the 
gammacoronavirus group too, which suggested a low 
prevalence in wild aquatic birds in USA, corroborating 
with our data [21]. As we know, migratory routes in 
Brazil are linked with flyways coming from North 
America and can confirm our findings [16]. Although 
our data have shown low prevalence when compared with 
studies before in Scandinavian Waterfowl, our data is 
clear when compared with North American data. In addi- 
tion, it is the first time that samples related to Gamma and 
deltacoronavirus group in free wild birds were reported in 
Brazil. 

The overall prevalence of CoV found in this study was 
0.8% (6/746 birds), with 3.7% (3/81) from Ibirapuera Park 
and 1.6% (3/192) from migratory birds captured in Lagoa do 
Peixe National Park. Other studies report an overall prevalence 
of CoV in aquatic and migratory birds between 0.33 and 18.5% 
[8, 17, 21, 30, 31]. This data reinforces that CoVs are circulat- 
ing across major migratory routes worldwide. 

Our report points to the circulation of coronaviruses 
among migratory and resident birds sampled from distinct 
regions of Brazil. Coronavirus RNA were detected in birds 
species of Anatidae family (Chinese goose) in Ibirapuera 
Park, Sao Paulo, southeast Brazil. Anatidae birds were al- 
ready described as hosts of coronavirus all over the globe. 
However, interestingly, we report for the first time the detec- 
tion in two sandpiper species (Calidris alba and Calidris 
fuscicollis) from Lagoa do Peixe National Park, a wintering 
area in southern Brazil. Although our data indicate a low 
prevalence of coronavirus in birds sampled, the genetic sim- 
ilarity between coronavirus detected at birds from distinct 
continents may indicate that wild birds are possibly carrying 
coronavirus throughout migratory routes as observed for avi- 
an influenza virus. We report the detection of samples related 
to delta- and gammacoronaviruses in different species of 
birds, some sampled in a park located inside the largest city 
in South America, indicating the circulation of coronavirus in 
urban areas (Table 3). In addition, all birds tested in this study 
were in apparently healthy populations. Further studies are 
required to determine the species and full genome similarity 
with other avian CoVs. However, our findings indicate that 
CoVs may be circulating among one of the major avian mi- 
gratory routes in the world. 
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